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ABSTRACT 

The reliability of quiet Sun magnetic field diagnostics based on the Fe I lines at 6302 A has 
been questioned by recent work. We present here the results of a thorough study of high-resolution 
multi-line observations taken with the new spectro-polarimeter SPINOR, comprising the 5250 and 
6302 A spectral domains. The observations were analyzed using several inversion algorithms, 
including Milne-Eddington, LTE with 1 and 2 components, and MISMA codes. We find that 
the line-ratio technique applied to the 5250 A lines is not sufficiently reliable to provide a direct 
magnetic diagnostic in the presence of thermal fluctuations and variable line broadening. In gen- 
eral, one needs to resort to inversion algorithms, ideally with realistic magneto-hydrodynamical 
constrains. When this is done, the 5250 A lines do not seem to provide any significant advantage 
over those at 6302 A. In fact, our results point towards a better performance with the latter 
(in the presence of turbulent line broadening). In any case, for very weak flux concentrations, 
neither spectral region alone provides sufficient constraints to fully disentangle the intrinsic field 
strengths. Instead, we advocate for a combined analysis of both spectral ranges, which yields a 
better determination of the quiet Sun magnetic properties. Finally we propose the use of two 
other Fe I lines (at 4122 and 9000 A) with identical line opacities that seem to work much better 
than the others. 

Subject headings: line: profiles - Sun: atmosphere - Sun: magnetic fields - Sun: photosphere 



1. Introduction 

The empirical investigation of quiet Surjf] mag- 
netism is a very important but extremely chal- 
lenging problem. A large (probably dominant) 
fraction of the solar magnetic flux resides in 
magnetic accumulations outside active regions, 
forming network and inter- network patch es (e.g., 
Socas-Navarro fe Sanchez Almeidal 120021 ). It is 



1 The National Center for Atmospheric Research 
(NCAR) is sponsored by the National Science Foundation. 
2 In this work, we use the term "quiet Sun" to refer to the 
solar surface away from sunspots and active regions. 



difficult to obtain conclusive observations of these 
structures, mainly because of two reasons. First, 
the size of the magnetic concentrations is much 
smaller than the spatial resolution capability of 
modern spectro-polarimetric instrumentation. Es- 
timates obtained with inversion codes yield typical 
values for the filling factor of the resolution cle- 
ment between ~1% and 30% . The interpretation 
of the polarization signal becomes non-trivial in 
these conditions and one needs to make use of de- 
tailed inversion codes to infer the magnetic field in 
the atmosphere. Second, the observed signals are 
extremely weak (typically below ~1% of the av- 
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erage continuum intensity), demanding both high 
sensitivity and high resolution. Linear polariza- 
tion is rarely observed in visible lines, so one is 
usually left with Stokes I and V alone. 

IStenflol ( 1973f ) proposed to use the pair of Fe I 
lines at 5247 and 5250 A which have very simi- 
lar excitation potentials and oscillator strengths 
(and, therefore, very similar opacities) but dif- 
ferent Lande factors, to determine the intrinsic 
field strength directly from the Stokes V line ra- 
tio. That work led to the subsequent populariza- 
tion of this spectral region for further studies of 
unresolved solar magnetic structures. Later, the 
pair of Fe I lines at 6302 A became the primary 
target of the Adva nced Stokes Polarimeter (ASP, 



Elmore et al.lll992T ) , mainly due to their lower sen- 



sitivity to temperature fluctuations. The success 
of the ASP has contributed largely to the currently 
widespread use of the 6302 A lines by the solar 
community. 

Recent advances in infrared spectro-polarimetric 
instrumentation now permit the routine obser- 
vation of another very interesting pair of Fe I 
lines, namely those at 15648 and 15653 A (here- 
after, t he 1.56 fjm lines). Examples are the 
works of iLin fe Rimmeld (| 19991 ): iKhomenko et al.l 
([20031 ) . The large Lande factors of these lines, 



combined with their very long wavelengths, re- 
sult in an extraordinary Zccman sensitivity. Their 
Stokes V profiles exhibit patterns where the a- 
components are completely split for fields stronger 
than ^400 G at typical photospheric conditions. 
They also produce stronger linear polarization. 
On the downside, this spectral range is accesi- 
ble to very few polarimeters. Furthermore, the 
1.56 /jm lines are rather weak in comparison with 
the above-mentioned visible lines. 

Unfortunately, the picture revealed by the new 
infrared data often differs drastically from what 
was being inferred from the 6302 A observa- 
tions (e .g., Lin fe Rimmeleill999llKhomenko et al 



20031: ISocas-Navarro fe Sanchez Almeidal 120021: 



sensitive to the weaker fields of the distribution 
whereas the visible lines will provide information 
on the stronger fields (see also the discussion about 
polarimetric signal incr ease in the 1.56 lines 
with weakening fields in Sanchez Almeida fe Litesl 
2000 ). This conjecture has been t ested recently 



bv lDommguez Cerdena et al.l (|2006f ) who modeled 



simultaenous observations of visible and infrared 
lines using unresolved magnetic inhomogeneities. 
A rece nt paper describing num erical simula- 



tions by iMartmez Gonzalez et al. ( 20061 ) casts 



Socas- Navarro fe Litesll2004T Dommguez Cerdefia et al 



2003) , particularly in the inter- network. ISocas-Navarro fe 
(|2003 ) proposed that the discrepancy in the field 
strengths inferred from the visible and infrared 
lines may be explained by magnetic inhomo- 
geneities within the resolution element (typically 
1"). If multiple field strengths coexist in the ob- 
served pixel, then the infrared lines will be more 



some doubts on the results obtained using the 
6302 A lines. Our motivation for the present 
work is to resolve this issue by observing simul- 
taneously the quiet Sun at 5250 and 6302 A. We 
know that unresolved magnetic structure might 
result in different field determinations in the visi- 
ble and the infrared, but the lines analyzed in this 
work are close enough in wavelengths and Zeeman 
sensitivities that one would expect to obtain the 
same results for both spectral regions. 

2. Methodology 

Initially, our goal was to observe simultane- 
ously at 5250 and 6302 A because we expected 
the 5250 A lines to be a very robust indicator of 
intrinsic field strength, which we could then use 
to test under what conditions the 6302 A lines 
are also robust. Unfortunately, as we show below, 
it turns out that in most practical situations the 
5250 A lines are not more robust than the 6302 A 
pair. This left us without a generally valid refer- 
ence frame against which to test the 6302 A lines. 
We then decided to employ a different approach 
for our study, namely to analyze the uniqueness 
of the solution obtained when we invert the lines 
and how the solutions derived for both pairs of 
lines compare to each other. In doing so, there are 
some subtleties that need to be taken into consid- 
eration. 

An inversion technique necessarily resorts on a 
number of physical assumptions on the solar at- 

Jjuosylieru n hi which the lines are formed and the 
Sanchez Almeida. . m , , 

(m general polarized) radiative transfer. 1ms im- 
plies that the conclusions obtained from applying 
a particular inversion code to our data may be 
biased by the modeling implicit in the inversion. 
Therefore, a rigorous study requires the analysis 
of solutions from a wide variety of inversion pro- 
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cedures. Ideally, one would like to cover at least 
the most frequently employed algorithms. For our 
purposes here we have chosen four of the most 
popular codes, spanning a wide range of model 
complexity. They are described in some detail in 
section [3] below. 

When dealing with Stokes inversion codes, 
there are two very distinct problems that the user 
needs to be aware of. First, it may happen that 
multiple different solutions provide satisfactory 
fits to our observations. This problem is of a very 
fundamental nature. It is not a problem with 
the inversion algorithm but with the observables 
themselves. Simply put, they do not carry suf- 
ficient information to discriminate among those 
particular solutions. The only way around this 
problem is to either supply additional observables 
(e.g., more spectral lines) or to restrict the al- 
lowed range of solutions by incorporating sensi- 
ble constraints in the physical model employed 
by the inversion. A second problem arises when 
the solution obtained does not fit the observed 
profiles satisfactorily. This can happen because 
the physical constraints in the inversion code are 
too stringent (and thus no good solution exists 
within the allowed range of parameters), or sim- 
ply because the algorithm happened to stop at a 
secondary minimum. This latter problem is not 
essential because one can always discard "bad" 
solutions (i.e., those that result in a bad fit to 
the data) and simply try again with a different 
initialization. In this work we are interested in 
exploring the robustness of the observables in as 
much as they relate to the former problem, i.e. 
the underlying uniqueness of the solution. We do 
this by performing a large number of inversions 
with random initializations and analyzing the de- 
pendence of the solutions with the merit function 
X 2 (defined below). 

Ideally, one would like to see that for small val- 
ues of x 2 ) all the solutions are clustered around 
a central value with a small spread (the behavior 
for large x 2 is not very relevant for our purposes 
here) . This should happen regardless of the inver- 
sion code employed and the pair of lines analyzed. 
If that were the case we could conclude that our 
observables are truly robust. Otherwise, one needs 
to be careful when analyzing data corresponding 
to that particular scenario. 



3. Observations and Data Reduction 

The observations used in this work were ob- 
tained during an observing run in March 2006 
with the Spectro-Polarimeter for Infrared and Op- 
tical R egions (SPINOR, see [s ocas-Navarro et al 



2006bh . attached to the Dunn Solar Telescope 
(DST) at the Sacramento Peak Observatory 
(Sunspot, NM USA, operated by the National 
Solar Observatory) . SPINOR allows for the simul- 
taneous observation of multiple spectral domains 
with nearly optimal polarimetric efficiency over a 
broad range of wavelengths. 

The high-order ada ptive optics system of the 
DST (|Ren et all [2001 was employed for image 
stabilization and to correct for atmospheric tur- 
bulence. This allowed us to attain sub-arcsecond 
spatial resolution during some periods of good see- 
ing. 

The observing campaign was originally devised 
to obtain as much information as possible on 
the unresolved properties of the quiet Sun mag- 
netic fields. In addition to the 5250 and 6302 A 
domains discussed here, we also observed the 
Mn I line sensitive to hv perfine structure effects 
((Lopez Ariste et al.|[2002l) at 5537 A and the Fe I 
line pair at 1.56 /im. 

SPINOR was operated in a configuration with 
four different detectors which were available at the 
time of observations (see Table [IJ: The Rock- 
well TCM 8600 infrared camera, with a format 
of 1024x1024 pixels, was observing the 1.56 /jm 
region. The SARNOFF CAM1M100 of 1024x512 
pixels was used at 5250 A. Finally, the two dual TI 
TC245 cameras of 256x256 pixels (the original de- 
tectors of the Advanced Stokes Polarimeter) were 
set to observe at 5537 and 6302 A. Unfortunately, 
we encountered some issues during the reduction 
of the 5537 A and 1.56 fim data and it is un- 
clear at this point whether or very not they are 
usable. Therefore, in the remainder of this paper 
we shall focus on the analysis of the data taken at 
5250 and 6302 A. The spectral resolutions quoted 
in the table are estimated as the quadratic sum 
of the spectrograph slit size, diffraction limit and 
pixel sampling. 

In order to have good spectrograph efficiency at 
all four wavelengths simultaneously we employed 
the 308.57 line mm -1 grating (blaze angle 52°), at 
the expense of obtaining a relatively low dispersion 
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Table 1 

spinor detector configuration 



Camera 


Wavelength 


Spectral resolution 


Spectral sampling 


Usable range 


Field of view 




(A) 


(mA) 


(mA) 


(A) 


along the slit (") 


ROCKWELL 


15650 


190 


150 


150 


187 


SARNOFF 


5250 


53 


31 


15 


145 


TI TC245 


5537 


40 


29 


5 


95 


TI TC245 


6302 


47 


24 


6 


95 



and spectral resolution (see Tabled] for details). 

Standard Hatfield and bias correction proce- 
duress were applied to the images. Subsequent 
processing included the removal of spectrum cur- 
vature and the alignment of both polarized beams, 
using for this purpose a pair of hairlines inserted 
across the slit. Calibration operations were per- 
formed to determine the SPINOR polarimetric re- 
sponse matrix by means of calibration optics lo- 
cated at the telescope beam exit port. In this 
manner we can decontaminate the datasets from 
instrumental polarization introduced by the po- 
larimeter. Finally, it is also important to consider 
the contamination introduced by the telescope. To 
this aim we obtained telescope calibration data 
with an array of linear polarizers situated over the 
DST entrance window. By rotating these polariz- 
ers to different angles, it is possible to feed light 
in known polarization states into the system. A 
cross-dispersing prism was placed in front of one 
of the detectors, allowing us to obtain calibration 
data simultaneously across the entire visible spec- 
trum. Details on th e proced ure may be found in 
Socas-Navarro et all (|2006bl ). 

In this paper we focus on two scan opera- 
tions near disk center, one over a relatively large 
pore (at solar heliocentric coordinates longitude 
-25.40, latitude -3.68) and the other of a quiet re- 
gion (coordinates longitude 0.01°, latitude -7.14°). 
The pore map was taken with rather low spatial 
resolution (~1.5") but exhibits a large range of 
polarization signal amplitudes. The quiet map, 
on the other hand, has very good spatial res- 
olution (~0.6") but the polarization signals are 
much weaker. The noise level, measured as the 
standard deviation of the polarization signal in 
continuum regions, is approximately 7xl0 -4 and 



5x10 4 times the average quiet Sun continuum 
intensity at 5250 and 6302 A, respectively. 

We used several different inversion codes for 
the various tests presented here, namely: SIR 
(Stokes Inversion based on Respon se functions, 
iRuiz Cobo fe del Toro Iniestal Il992j ): MELANIE 
(Milne-Eddington Line Analysis using a Nu- 
merical Inversion Engine) and LILIA (LTE In- 
versio n based on the L orien Iterative Algo- 
rithm, [SocaslN^r3[200H); and MISMA (Micro- 
Struc tured Magnetic Atmosphere. ISanchez Almeida 
19971 ). The simplest of these algorithms is 



MELANIE, which implements a Milne-Eddington 
type o f inversion similar to that of Skumanich fe Litesl 
(|l987l ). The free parameters considered include a 
constant along the line of sight magnetic field vec- 
tor, magnetic filling factor, line-of-sight velocity 
and several spectral line parameters that repre- 
sent the thermal properties of the atmosphere 
(Doppler width AAc line-to-continuum opacity 
ratio ?7o, source function S and damping a). The 
Fe I lines at 6302 A belong to the same multiplct 
and their 770 are related by a constant factor. As- 
suming that the formation height is similar enough 
for both lines we can also consider that they have 
the same AA,d, S and a. In this manner the same 
set of free parameters can be used to fit both lines. 
In the case of the 5250 A lines we only invert the 
Fe I pair and assume that both lines have identical 
opacities tjq. 

SIR considers a model atmosphere defined by 
the depth stratification of variables such as tem- 
perature, pressure, magnetic field vector, line-of- 
sight velocity and microturbulence. Atomic pop- 
ulations are computed assuming LTE for the vari- 
ous lines involved, making it possible to fit obser- 
vations of lines from multiple chemical elements 
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with a single model atmosphere that is common 
to all of them. Unlike MELANIE, one can pro- 
duce line asymmetries by incorporating gradients 
with height of the velocity and other parameters. 
LILIA is a different implementation of the SIR al- 
gorithm. It works very similarly with some prac- 
tical differences that are not necessary to discuss 
here. 

MISMA is another LTE code but has the capa- 
bility to consider three atmospheric components 
(two magnetic and one non-magnetic) that are in- 
terlaced on spatial scales smaller than the photon 
mean free path. Perhaps the most interesting fea- 
ture of this code for our purposes here is that it 
implements a number of magneto-hydrodynamic 
(MHD) constrains, such as momentum, as well as 
mass and magnetic flux conservation. In this man- 
ner it is possible to derive the full vertical strat- 
ification of the model atmosphere from a limited 
number of free parameters (e.g., the magnetic field 
and the velocity at the base of the atmosphere). 

In all of the inversions presented here we em- 
ployed the same set of atomic line parameters, 
which are listed in Table [2] 

4. Results 

Figures [T] and [2] show continuum maps and re- 
constructed magnetograms of the pore and quiet 
Sun scans in the 5250 A region. Notice the much 
higher spatial resolution in the quiet Sun observa- 
tion (Fig[l|. Similar maps, but with a somewhat 
smaller field of views, can be produced at 6302 A. 

The first natural step in the analysis of these 
observations, before even considering any inver- 
sions, is to calculate the Stokes V amplitude ra- 
tio of the Fe I line 5250.2 to 5247.0 A (hereafter, 
the line ratio). One would expect to obtain a 
rough idea of the intrinsic magnetic field strength 
from this value alone. A simple calibration was 
derived by taking the thermal stratification of 
the Harvar d-Smithsonian Refere nce Atmosphere 
(HSRA, see lGingerich et al1ll971 ) and adding ran- 
dom (depth-independent) Gaussian temperature 
perturbations with an amplitude (standard devia- 
tion) of ±300 K, different magnetic field strengths 
and a fixed macroturbulence of 3 km s _1 (this 
value corresponds roughly to our spectral resolu- 
tion). Line-of-sight gradients of temperature, field 
strength and velocity are also included. Figure [3] 



shows the line ratio obtained from the synthetic 
profiles as a function of the magnetic field em- 
ployed to synthesize them. The line ratio is com- 
puted simply as the peak-to-peak amplitude ra- 
tio of the Stokes V profiles. Notice that when 
the same experiment is carried out with a variable 
macroturbulent velocity the scatter increases con- 
siderably, even for relatively small values of up to 
lkms -1 (right panel). The syntheses of Fig[3]con- 
sider the partial blends of all 6 lines in the 5250 A 
spectral range. 

Figure 2] shows the observed line ratios for 
both scans. According to our calibration (see 
above), ratios close to 1 indicate strong fields of 
nearly (at least) ~2 kG, whereas larger values 
would suggest the presence of weaker fields, down 
to the weak-field saturation regime at (at most) 
~500 G corresponding to a ratio of 1.5. Fig- 
ure S] is somewhat disconcerting at first sight. 
The pore exhibits the expected behavior with 
strong ~2 kG fields at the center that decrease 
gradually towards the outer boundaries until it 
becomes weak. The network and plage patches, 
on the other hand, contain relatively large areas 
with high ratios of 1.4 and even 1.6 at some lo- 
cations. This is in sharp contrast with the strong 
fields (~1.5 kG) that one w ould expect in net- 
work and plage re gions (e.g.. IStenflo et al.lll984t 



Sanchez Almeida & Lite sboOoilBellot Rubio et al 
2000). 

In view of these results we carried out inversions 
of the Stokes / and V profiles of the spectral lines 
in the 5250 A region emergent from the pore. We 
used the code LILIA considering a single magnetic 
component with a constant magnetic field embed- 
ded in a (fixed) non-magnetic background, taken 
to be the average quiet Sun. The resulting field- 
ratio scatter plot is shown in Figure O Note that 
the scatter in this case is much larger than that 
obtained with the HSRA calibration above. It is 
important to point out that the line ratio depicted 
in the figure is that measured on the synthetic 
profiles. Therefore, the scatter cannot be ascribed 
to inaccuracies of the inversion. As a verification 
test we picked one of the models with kG fields 
that produced a line ratio of ~ 1.4 and synthesized 
the emergent profiles with a different code (SIR), 
obtaining the same ratio. We therefore conclude 
that, when the field is strong, many different at- 
mospheres are able to produce similar line ratios 
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Table 2 
Atomic line data 





Wavelength 


Excit. Potential 


log(ff/) 


Lower level 


Upper level 


Ion 


(A) 


(eV) 




configuration 


configuration 


Fe I 


4121.8020 


2.832 


-1.300 


3 P 2 


3 F 3 


Cr II 


5246.7680 


3.714 


-2.630 


4 Pl/2 


4 P 3/2 


Fe I 


5247.0504 


0.087 


-4.946 


5 D 2 


7 D 3 


Til 


5247.2870 


2.103 


-0.927 


5 F 3 


5 F 2 


CrI 


5247.5660 


0.961 


-1.640 


5 D 


5 Pi 


Fe I 


5250.2089 


0.121 


-4.938 


5 D 




Fe I 


5250.6460 


2.198 


-2.181 


5 P 2 


5 P 3 


Fe I 


6301.5012 


3.654 


-0.718 


5 P 2 


5 D 2 


Fe I 


6302.4916 


3.686 


-1.235 


5 Pi 


5 D 


Fe I 


8999.5600 


2.832 


-1.300 


3 P 2 


3 P 2 




10 20 30 40 10 20 30 40 

Arc-seconds Arc-seconds 



Fig. 1. Continuum intensity map (left) and reconstructed magnetogram (right) of the quiet Sun scan near 
5250 A. The magnetogram shows the integrated Stokes V signal over a narrow bandwidth on the red lobe of 
the Fe I 5250.2 A line. In both panels, values are referred to the average quiet Sun intensity. The angular 
scale on the vertical axis measures the distance from the upper hairline. 
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Fig. 2. — Continuum intensity map (left) and reconstructed magnetogram (right) of the pore scan near 
5250 A. The magnetogram shows the integrated Stokes V signal over a narrow bandwidth on the red lobe of 
the Fe I 5250.2 A line. In both panels, values are referred to the average quiet Sun intensity. The angular 
scale on the vertical axis measures the distance from the upper hairline. 
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Fig. 3. Line ratio calibration using synthetic Stokes profiles of Fe I 5250 A emergent from the HSRA model 
atmosphere after adding random temperature perturbations and magnetic fields. Left: Macroturbulence was 
kept fixed at 3 km s _1 . Right: Random macroturbulence varies between and 1 km s _1 . 
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Fig. 4. — Amplitude ratio of Stokes V in Fe I 5250.2 to Fe I 5247.0 A. Left: Network patches observed in 
the quiet Sun scan. Right: A small pore (centered on coordinates [10,-30] approximately) and surrounding 
plage. Ratios close to 1 are indicative of a strong field of nearly 2 kG. Ratios close to 1.6 correspond to the 
weak-field (up to ~500 G) regime. Black areas exhibit a circular polarization amplitude smaller than 1% 
and have not been inverted. 



if realistic thermodynamic fluctuations and tur- 
bulence are allowed in the model. For weak fields, 
the ratio tends to a value of ~1.5 without much 
fluctuations. 

Accepting then that we could no longer rely on 
the line ratio of Fe I 5247.0 and 5250.2 A as an 
independent reference to verify the magnetic fields 
obtained with the 6302 A lines, we considered the 
result of inverting each spectral region separately. 
Figure [6] depicts the scatter plot obtained. At the 
center of the pore, where we have the strongest 
fields (right-hand side of the figure), there is a very 
good correlation between the results of both mea- 
surements. However, those points lay systemati- 
cally below the diagonal of the plot. This may be 
explained by the different "formation heights" of 
the lines. The Fe I lines at 5250 A generally form 
somewhat higher than those at 6302 A. If the field 
strength decreases with height, one would expect 
to retrieve a slightly lower field strength when us- 
ing 5250 A. Unfortunately, the correlation breaks 
down for the weaker fields. In Figure [7] we can 
see that both sets of lines yield approximately the 
same field strengths (with a slightly lower values 
for the 5250 inversions, as discussed above) for lon- 
gitudinal fluxes above some ^500 G. Below this 
limit our diagnosis is probably not sufficiently ro- 
bust. 



1.6 




1.0 I , , , , I , , , , I , , , , I , , , , 

500 1000 1500 2000 

Field strength 



Fig. 5. — Line ratio calibration using synthetic 
Stokes profiles obtained in the inversion of the 
pore. 
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1DO0 2000 
Field (G) from 6302 inversions 



Fig. 6. — Intrinsic field strengths determined from 
inverting the lines at 5250 (ordinates) and 6302 
(abscissa) observed in the pore. The solid line 
shows the diagonal of the plot. 




If instead of the intrinsic field we consider the 
longitudinal magnetic flux, we obtain a fairly good 
agreement between both spectral regions. The 
Milnc-Eddington inversions with MELANIE yield 
a Pearsons correlation coefficient of 0.89. In prin- 
ciple, the agreement is somewhat worse for the 
LILIA inversions, with a correlation coefficient of 
0.60. However, we have found that this is due to 
a few outlayer points. Removing them results in 
a correlation coefficient of 0.82. 

The situation becomes more complicated in the 
network. The results of inverting a network patch 
with MELANIE and LILIA can be seen in Fig- 
ure[8l The inversions with MELANIE (upper pan- 
els) do not agree very well with each other. The 
correlation coefficient is only 0.23. The 5250 A 
map (upper right panel in the figure) looks consid- 
erably more noisy and rather homogeneous, com- 
pared to the 6302 A map (upper left) . The LILIA 
inversions (lower panels) exhibit somewhat bet- 
ter agreement (correlation is 0.65), but again look 
very noisy at 5250 A. 

Figure [9] depicts average profiles over a net- 
work patch. The 5250 line ratio for this profile 
is 1.21. We started by exploring the uniqueness 
of the magnetic field strength inferred with the 
simplest algorithm, MELANIE. Each one of these 
average profiles was inverted 100 times with ran- 
dom initializations. The results are presented in 
Figures [inland [Til which show the values obtained 
versus the goodness of the fit, defined as the merit 




1000 

Longitudinal flux (G) 









b J 















Fig. 7. — Differences in intrinsic field strengths 
measured at 6302 and 5250 A as a function of the 
average longitudinal magnetic flux density. The 
solid line represents the median value over 100 G 
bins. The dashed lines represent the standard de- 
viation of the points in each bin. 



Fig. 8. — Longitudinal flux density from the inver- 
sions of a network patch at 6302 (left) and 5250 A 
(right). Upper panels: Inversions with the Milnc- 
Eddington code MELANIE. Lower panels: Inver- 
sions with the LTE code LILIA. Spatial scales are 
in pixels. 
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function x 2 , which in this work is defined as: 



1 N p 



obs 



r syn\2 



(i) 



where N p is the number of wavelengths and <7; 
have been taken to be 10~ 3 , so that a value of 
X 2 = l would represent on average a good fit at the 
10 -3 level. The average profiles inverted here have 
a much lower noise (near 10 ) and thus it is some 
times possible to obtain x 2 smaller than 1. 

The x 2 represented in the plots is the one cor- 
responding to the Stokes V profile only (although 
the inversion codes consider both / and V, but I is 
consistently well reproduced and does not help to 
discriminate among the different solutions). Most 
of the fits correspond to kG field, indicating that 
inversions of network profiles are very likely to 
yield high field strengths. However, there exists 
a very large spread of field strength values that 
provide reasonably good fits to the observed data. 
This is especially true for the 5250 A lines, for 
which it is possible to fit the observations vir- 
tually equally well with fields either weaker than 
500 G or stronger than 1 kG. In the case of 6302 A 
the best solutions are packed around ~1.5 kG, al- 
though other solutions of a few hecto-Gauss (hG) 
are only slightly worse than the best fit. 

It could be argued that Milnc-Eddington in- 
versions are too simplistic to deal with network 
profiles, since they are known to exhibit fairly 
strong asymmetries (both in area and in ampli- 
tude) that cannot be reproduced by a Milne- 
Eddington model. With this consideration in 
mind, we made a similar experiment using the 
LILIA and SIR codes. Figures [12] and [13] show 
the results obtained with LILIA. The magnetic 
and non-magnetic atmospheres have been forced 
to have the same thermodynamics in order to re- 
duce the possible degrees of freedom (on the down- 
side, this introduces an implicit assumption on 
the solar atmosphere). Again, the 6302 A lines 
seem to yield somewhat more robust inversions. 
The best fits correspond to field strengths of ap- 
proximately 1.5 kG, with weaker fields delivering 
somewhat lower fit quality. The 5250 A lines give 
nearly random results (although they tend to be 
clustered between 500 and 800 G there is a tail 
of good fits with up to almost 1400 G). Similar 
results are obtained using SIR. In order to give 
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A - 6302 (A) 

Fig. 9. — Average observed Stokes V profiles in 
a network patch (located approximately around 
coordinates [20,7] in Figure [T]). Ordinate values 
are related to the average quiet Sun continuum 
intensity. The 5250 line ratio (the relevant Fe I 
lines are marked with arrows) is 1.21. 




Fig. 10. — Representation of the solutions from 
100 different inversions of the 5250 A region with 
random initializations obtained using the Milne- 
Eddington code MELANIE. The field strength in- 
ferred by the inversion is plotted versus the quality 
of the fit, as measured by the merit function x 2 ■ 
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an idea of what the different \ 2 values mean, we 
present some of the fits in Fig IT4l 

The thermal stratifications obtained in the 
6302 A inversions are relatively similar, although 
the weaker fields require a hotter upper photo- 
sphere than the stronger fields (see Fig [T5|> . On 
the other hand, the 5250 A inversions do not 
exhibit a clear correlation between the magnetic 
fields and temperature inferred fFigfTo]). 

The M1SMA code was also able to find good 
solutions with either weak or strong fields. The 
smallest \ 2 values correspond systematically to 
kG fields, but there are also some reasonably good 
fits (x 2 — 1) obtained with weak fields of ^500 G 
(Figs [T7] and US]). However, we found that all the 
weak-field solutions for 6302 A have a tempera- 
ture that increases outwards in the upper photo- 
sphere (Fig [T7|) . This might be useful to discrim- 
inate between the various solutions. The 5250 A 
lines, on the other hand, do not exhibit this be- 
havior (Fig ITS)) . 

In principle it would seem plausible to discard 
the models with outward increasing temperature 
using physical arguments. This would make us 
conclude from the 6302 A lines that the fields are 
actually very strong, between ~1.5 and 2.5 kG (it 
would not be possible to draw similar conclusions 
from the 5250 A lines). In any case, it would be 
desirable to have a less model-dependent measure- 
ment that could be trusted regardless of what the 
thermal stratification is. 

Interestingly, when we invert all four Fe I lines 
simultaneously, both at 5250 and 6302 A, the 
weak-field solutions disappear from the low-x 2 re- 
gion of the plot and the best solutions gather be- 
tween approximately 1 and 1.4 kG (see Fig [19]). 

Socas-Navarro et all ( 2006a[ ) present a list of 



spectral lines with identical excitation potentials 
and oscillator strengths. We decided to test one 
of the most promising pairs, namely the Fe I lines 
at 4122 and 9000 A. The choice was made based 
on their equivalent widths, Lande factors and also 
being reasonably free of blends in the quiet solar 
spectrum. Simulations similar to those in Figure [3] 
showed that the line-ratio technique is still inca- 
pable of retrieving an unambiguous field strength 
due to the presence of line broadening. Also, even 
if the line opacities are the same, the continuum 
opacities are significantly different at such dis- 



Fig. 11. — Representation of the solutions from 
100 different inversions of the 6302 A region with 
random initializations obtained using the Milnc- 
Eddington code MELANIE. The field strength in- 
ferred by the inversion is plotted versus the quality 
of the fit, as measured by the merit function x 2 . 



Fig. 12. — Representation of the solutions from 
100 different inversions of the 5250 A region with 
random initializations obtained using the LTE 
code LILIA. The field strength inferred by the in- 
version is plotted versus the quality of the fit, as 
measured by the merit function x 2 . 
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Fig. 14. — Fits obtained with MELANIE (left) and LILIA (right) to the average network profile observed at 
5250 (upper four panels) and 6302 A (lower four panels). For each case we show the best fit (smaller x 2 ) 
and a not so good one. 
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Fig. 13. — Representation of the solutions from 
100 different inversions of the 6302 A region with 
random initializations obtained using the LTE 
code LILIA. The field strength inferred by the in- 
version is plotted versus the quality of the fit, as 
measured by the merit function x 2 . 



Fig. 16. — Temperature stratification of the mod- 
els obtained from the 5250 A inversions with 
LILIA. Solid line: Models that include kG fields. 
Dashed line: Models with sub-kG field. 
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Fig. 15. — Temperature stratification of the mod- 
els obtained from the 6302 A inversions with 
LILIA. Solid line: Models that include kG fields. 
Dashed line: Models with sub-kG field. 



Fig. 17. — Inversions of the 6302 A region. Aver- 
age magnetic field in the lower photosphere of the 
MISMA component that harbors more flux, as a 
function of the quality of the fit x 2 - Asterisks: So- 
lutions where the temperature decreases outwards 
in the upper photosphere. Diamonds: Solutions 
where the temperature increases outwards in the 
upper photosphere. 
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Fig. 18. — Inversions of the 5250 A region. Aver- 
age magnetic field in the lower photosphere of the 
MISMA component that harbors more flux, as a 
function of the quality of the fit \ 2 ■ Asterisks: So- 
lutions where the temperature decreases outwards 
in the upper photosphere. Diamonds: Solutions 
where the temperature increases outwards in the 
upper photosphere. 




Fig. 19. — Inversions of combined 5250 and 6302 A 
profiles. Average magnetic field in the lower pho- 
tosphere of the MISMA component that harbors 
more flux, as a function of the quality of the fit 

x 2 - 



parate wavelengths, but in any case the line opac- 
ity is much stronger than the background contin- 
uum where the Stokes V lobes are formed. 

When we tested the robustness of inversion 
codes applied to this pair of lines, we obtained 
extremely reliable results as described below. Un- 
fortunately we do not have observations of these 
two lines and therefore resorted on synthetic pro- 
files to perform the tests. We used a 2-component 
reference model, where both components have the 
thermal structure of HSRA. The magnetic como- 
nent has an arbitrary magnetic field and line-of- 
sight velocity. The field has a linear gradient and 
goes from ~1.6 kG at the base of the photosphere 
to 1 kG at continuum optical depth log(r5ooo)=- 
4. The velocity field goes from 2.5 km s _1 to 
0.4 km s _1 at log(r5ooo)=-4. The magnetic com- 
ponent has a filling factor of 0.2. The reference 
macroturbulence is 3 km s _1 , which is roughly the 
value inferred from our observations. The spec- 
tra produced by this model at 5250 and 6302 A 
are very similar to typical network profiles. We 
considered the reference 4122 and 9000 A profiles 
as simulated observations and inverted them with 
two different methods. 

In order to make the test as realistic as possible, 
we gave the inversion code a somewhat erroneous 
non-magnetic profile. Specifically, we multiplied 
the profile from the reference non-magnetic com- 
ponent by a factor of 0.9 and shifted it 1 pixel to- 
wards the red. The wavelength shift corresponds 
to roughly 0.36 km s _1 in both regions. This dis- 
torted non-magnetic profile was given as input to 
the inversion codes. We inverted the reference pro- 
files using LILIA with 100 different initializations. 
Only ^30% of the inversions converged to a rea- 
sonably low value of \ 2 , with the results plotted in 
Figure [201 We can see that the inversions are ex- 
tremely consistent over a range of x 2 much larger 
than in the previous cases. In fact, none of the so- 
lutions are compatible with weak fields, suggesting 
that these lines are much better at discriminating 
intrinsic field strengths. 

A much more demanding verification for the di- 
agnostic potential of these new lines is to use a sim- 
pler scenario in the inversion than in the synthesis 
of the reference profiles, incorporating typical un- 
certainties in the calculation. After all, the real 
Sun will always be more complex than our simpli- 
fied physical models. Thus, inverting the reference 
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profiles with the Milne-Eddington code is an ap- 
propriate test. For this experiment we not only 
supplied the same "distorted" non-magnetic pro- 
file as above, but we also introduced a systematic 
error in the log(<?/) of the lines. We forced the 
opacity of the 9000 A line to be 20% lower than 
that of 4122 (instead of taking them to be iden- 
tical, as their tabulated values would indicate). 
Again, we performed 100 different inversions of 
the reference set of profiles with random initial- 
izations. In this case the inversion results are as- 
tonishingly stable, with 98 out of the 100 inver- 
sions converging to a x 2 within 15% of the best 
fit. The single- valued magnetic field obtained for 
those 98 inversions has a median of 1780 G with a 
standard deviation of only 3 G. The small scatter 
does not reflect the systematic errors introduced 
by several factors, including: a)the inability of the 
Milne-Eddington model to reproduce the compar- 
atively more complex referece profiles; b)the arti- 
ficial error introduced in the atomic parameters of 
the 9000 A line; c)the distortion (scale and shift) 
of the non-magnetic profile provided to the inver- 
sion code. 

A final caveat with this n ew pair is that, even 



though the synthetic atlas of iSocas-Navarro et al 



(|2006al ) indicates that the 9000 A line Stokes V 
profile is relatively free of blends, this still needs 
to be confirmed by observations (there is a very 
prominent line nearby that may complicate the 
analysis otherwise). 

5. Conclusions 

The ratio of Stokes V amplitudes at 5250 and 
5247 A is a very good indicator of the intrinsic 
field strength in the absence of line broadening, 
e.g. due to turbulence. However, line broadening 
tends to smear out spectral features and reduce 
the Stokes V amplitudes. This reduction is not 
the same for both lines, depending on the profile 
shape. If the broadening could somehow be held 
constant, one would obtain a line-ratio calibration 
with very low scatter. However, if the broadening 
is allowed to fluctuate, even with amplitudes as 
small as 1 km s _1 , the scatter becomes very large. 
Fluctuations in the thermal conditions of the at- 
mosphere further complicate the analysis. This 
paper is not intended to question the historical 
merits of the line-ratio technique, which led re- 



searchers to learn that fields seen in the quiet Sun 
at low spatial resolution are mostly of kG strength 
with small filling factors. However, it is impor- 
tant to know its limitations. Otherwise, the inter- 
pretation of data such as those in Figure 2] could 
be misleading. Before this work, most of the au- 
thors were under the impression that measuring 
the line ratio of the 5250 A lines would always 
provide an accurate determination of the intrinsic 
field strength. 

With very high-resolution observations, such as 
those expected from t he Advanced Tec hnology So- 
lar Telescope (ATST, iKeil et alJliool or the Hin- 
ode satellite, there is some hope that most of the 
turbulent velocity fields may be resolved. In that 
case, the turbulent broadening would be negligi- 
ble and the line-ratio technique would be more 
robust. However, even with the highest possible 
spatial resolution, velocity and temperature fluc- 
tuations along the line of sight will still produce 
turbulent broadening. 

From the study presented here we conclude 
that, away from active region flux concentrations, 
it is not straightforward to measure intrinsic field 
strengths from either 5250 or 6302 A observa- 
tions taken separately. Weak-flux internetwork 
observations would be even more challenging, 
as dem onstrated recently by iMartinez Gonzalez 
( 2007 ). Surprisingly enough, the 6302 A pair of 
Fe I lines is more robust than the 5250 A lines in 
the sense that it is indeed possible to discriminate 
between weak and strong field solutions if one is 
able to rule out a thermal stratification with tem- 
peratures that increase outwards. Even so, this is 
only possible when one employs an inversion code 
that has sufficient MHD constrains (an example is 
the MISMA implementation used here) to reduce 
the space of possible solutions. 

The longitudinal flux density obtained from in- 
versions of the 6302 A lines is better determined 
than those obtained with 5250 A. This happens 
regardless of the inversion method employed, al- 
though using a code like LILIA provides better 
results than a simpler one such as MELANIE. The 
best fits to average network profiles correspond to 
strong kG fields, as one would expect. 

An interesting conclusion of this study is that 
it is possible to obtain reliable results by invert- 
ing simultaneous observations at both 5250 and 
6302 A. Obviously this would be possible with rel- 
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atively sophisticated algorithms (e.g., LTE inver- 
sions) but not with simple Milne-Eddington inver- 
sions. 

The combination of two other Fe I lines, namely 
those at 4122 and 9000 A, seems to provide a much 
more robust determination of the quiet Sun mag- 
netic fields. Unfortunately, these lines are very dis- 
tant in wavelength and few spectro-polarimeters 
are capable of observing them simultaneously. Ex- 
amples of instrument with this capability are the 
currently operational SPINOR and THEMIS, as 
well as the planned ATST and GREGOR. De- 
pending on the evolution time scales of the struc- 
tures analyzed it may be possible for some other 
instruments to observe the blue and red lines al- 
ternatively. 

This work has been partially funded by the 
Spanish Ministerio de Education y Ciencia through 
project AYA2004-05792 
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Fig. 20. — Inversions of the proposed Fe I 4122 
and 9000 A profiles. Average magnetic field in 
the lower photosphere of the magnetic component 
as a function of the quality of the fit \ 2 ■ 



